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Abstract . The production of a single neutral Higgs boson h through (loop-induced) 
77 collisions is explored in the context of the linear colliders within the general Two- 
Higgs-Doublet Model (2HDM). Two different mechanisms are analyzed: on the one 
hand, the scattering 77 ^ /i of two real photons in a 77 collider; on the other, the more 
traditional mechanism of virtual photon fusion, e^e~ — > e^e~7*7* — > e^e~ -|- h. Owing 
to the peculiar properties of the Higgs boson self-interactions within the general 2HDM, 
we find that the overall production rates can be boosted up significantly, provided 
the charged Higgs mass is not too heavy. For example, if Mjj± > 100 GeV and, in 
addition, M^o falls in the ballpark of the LEP bound on the SM Higgs mass up to a 
few hundred GeV, the cross-sections may typically render (a-y^^/i) ~ 0.1 — Ipb and 
cr(e+e- e+e-hO) < 0.01 pb - in both cases well above the SM prediction. Although 
for -Mjj± > 300 GeF the rates become virtually insensitive to the Higgs boson self- 
couplings, a significant tail of non-SM effects produced by the combined contribution 
of the Yukawa couplings and gauge bosons could still reveal a smoking gun. 
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1 Introduction 



The Higgs mechanism is the most fundamental hngering issue that remains experimentally unset- 
tled in Particle Physics. It is difficult to overemphasize that this issue stands right in the core of our 
present understanding of the Standard Model (SM). However, we cannot exclude that the Higgs 
sector is larger than expected, the most paradigmatic extension being the Minimal Super symmetric 
Standard Model (MSSM) [1], which involves two doublets of complex scalar fields. The physical 
spectrum consists of two charged states, H^, two neutral CP-even states hP,H^ (with masses 
M^o < Mfjo) and one CP-odd state [2]. Let us recall that the self-interactions of the SUSY 
Higgs bosons are rather inconspicuous, in the sense that they cannot be enhanced as compared 
to the ordinary gauge interactions and, therefore, do not present a very distinctive phenomenol- 
ogy. The bulk of the enhancing capabilities of the MSSM Lagrangian resides, instead, in the rich 
structure of Yukawa couplings between Higgs bosons and quarks or between quarks, squarks and 
chargino-neutralinos. The stupendous phenomenological opportunities associated to these super- 
symmetric structures are well-known since long ago (cf. [3,4]) and have been continuously updated 
in the literature (for recent reviews, see e.g. [5]). 

On the other hand, we should be prepared for alternative forms of Higgs boson physics of a 
more generic kind, whose potential implications can be equally outstanding and nevertheless be 
concentrated on very different sectors of the model. This could e.g. be the case of the general 
(unconstrained) Two-Higgs-Doublet Model (2HDM), where again two doublets of complex scalar 
fields are introduced, leading to a similar physical spectrum h^,H^,A^,H^, but without being 
subdued by the severe restrictions enforced by the supersymmetric transformations. The result 
is a Higgs potential with a collection of Higgs boson self-interactions which, in contrast to the 
SUSY case, can be highly enhanced in comparison to the gauge couplings. We refer the reader to 
Ref. [2,5] for further details. 

Let us assume that the LHC unveils a neutral Higgs boson. An essential part of the very process 
of identification will be to disclose just whether such particle is actually the neutral SM Higgs boson, 
a neutral member of a SUSY extension of the SM (typically the MSSM), or rather a generic neutral 
Higgs boson of a non-SUSY alternative setup, such as e.g. the general 2HDM. In this task, the 
complementary help of the future linear e~^e~ colliders (linac) [6], such as the ILC and CLIC, 
can play a momentous role to unravel the ultimate nature of the purported Higgs boson scalar 
particle(s) presumably produced at the LHC. There are many studies in the literature supporting 
this fact. For instance, the trilinear (3H) couplings have been investigated phenomenologically 
in TeV-class linear colliders in [7-10] through the double-Higgs strahlung process e''~e~ — > HHZ 
or the WW double-Higgs fusion mechanism e'''e~ — > H''"H~z4l7e. Unfortunately, the cross-section 
turns out to be rather small both in the SM and in the MSSM. Quite in contrast, it has recently 
been shown that the general 2HDM can provide cross-sections two to three orders of magnitude 
larger within the same experimental setup [11,12]. 

Closely connected to the physics of the linear colliders will be the physics of the 77 colliders [13] . 
As is well-known, a collider of this sort can be optionally realized from a linac by the process of 
backward Compton scattering between laser photons and the linac leptons. Not surprisingly, a 
clean machine as a 77 collider should enable us to probe the most sensitive theoretical structures 
of gauge theories, and certainly the Higgs sector is a most preeminent one. 

Currently, a renewed thrust of theoretical activity has been invested in double Higgs production 
in 77 collisions within the general 2HDM [14] - see also [15,16] for earlier related work, and [17,18] 
for single and double Higgs production in the SM and the MSSM. In this Letter, we wish to further 
explore the Higgs boson self-interactions in the general 2HDM by focusing on the process of single 
neutral Higgs boson production 77 — > /i = h*^, H*^, A*^ in the context of both 77 real scattering and 
77 virtual fusion in e"*" e~ colliders (see Fig.[T]). 
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2 Loop-induced 77/i interactions within the 2HDM: general fea- 
tures and computational setup 

Let us recall that the general 2HDM [2] is obtained by canonically extending the SM Higgs sector 
with a second SUi{2) doublet carrying weak hypercharge Y = +1, so that it contains 4 complex 
scalar fields. The free parameters in the most general, CP-conserving, 2HDM potential can be 
expressed in terms of the masses of the physical Higgs particles, M^o, Mjjo, Mpp, Mjj±, the ratio 
tan/3 = V2/V1 of the two VEV's giving masses to the up- and down-like quarks, the mixing angle 
a between the two CP-even states, and, finally, the coupling A5 which cannot be absorbed in any 
of the previous quantities^. In turn, the possible 2HDM coupling patterns in the Higgs-fermion 
sector are commonly sorted out as follows: i) type-I models, in which only one Higgs doublet 
couples to fermions, whereas the other doublet does not; and ii) type-H models, wherein a doublet 
couples only to down-like fermions and the other doublet only to up-like fermions. In either way 
one may avoid the appearance of dangerous (tree-level) Flavor Changing Neutral Current (FCNC) 
processes [2]. The MSSM Higgs sector is actually a type-II one, but of a very restricted sort 
(enforced by SUSY invariance) [1]. 

On top of that a number of important restrictions, emerging from either the available experi- 
mental data and the theoretical consistency of the model, must be taken into account in order to 
obtain a more realistic output. Although we have already described these constraints in [11, 12] 
(see also [19]), we will introduce some qualifications here. To start with, there are (additional) 
stringent constraints coming from (one-loop induced) low-energy FCNC processes, mainly from 
the charged Higgs boson contributions to B{h 57) [20], which require Mjj± > 295 GeV (for 
tan/3 ^ 1) in type-II models. Let us emphasize that this bound does not apply to type-I models 
since for them the charged Higgs couplings to fermions are proportional to cot (5 and hence the 
loop contributions are highly suppressed at large tan (3. Furthermore, the approximate SU{2) cus- 
todial symmetry severely restricts the radiative corrections to the p parameter from the 2HDM 
degrees of freedom; experimentally \6 p\ < 10~^ [21]. Moreover, there are of course the bounds 
stemming from the unsuccessful Higgs boson searches at LEP and the Tevatron [21]. Besides, a 
very important set of conditions emerges from the unitarity constraints. A substantial number 
of studies are devoted to this subject in the literature [22,23], although their conclusions are not 
always fully coincident. Alternatively, one can stick to a less restrictive (albeit well-motivated) 
framework based on requiring that none of the triple and quartic Higgs boson self-couplings (in 
the mass-eigenstate basis) should be larger than the upper value of the corresponding couplings 
in the SM. Ultimately, this condition is grounded on the Lee-Quigg-Thacker unitarity bound [24] 
on the SM Higgs boson mass. In this way, we obtain a double set of conditions that significantly 
harness the size of the 3H and 4H couplings in the 2HDM: 
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(— e being the electron charge and Oyy the weak mixing angle). In the following, we will discuss our 
results by taking into account the conditions ([I]) and ([2]), and we will briefly compare them with 
the restrictions derived in Ref. [22]. Furthermore, we shall impose that the EW vacuum is stable. 



^Throughout the paper, we use the notation and conventions of Ref. [11], to which we refer the reader for further 
details. Here, in contrast to that reference, we leave As as a fully independent parameter. 
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(a) 



(b) 



Figure 1: Generic Feynman diagrams describing the single Higgs production process through the mechanisms of 
a) direct scattering of a real photon pair and b) virtual 77 fusion in e'''e~ collisions. The grey blobs stand for the 
generic loop-induced 77/1 interaction at any order in perturbation theory. 



which is tantamount to say that we demand the quartic interaction terms in the potential not 
to give negative contributions producing an unbounded potential from below [25]. This condition 
leads to 

Ai > 0, A2 > 0, 
VA1A2 + A3 + Min(0, A4 + A5, A4 - A5) > 0, (3) 

where the parameters Aj are defined in terms of the Aj ones [11] as follows: 

Ai = 2(Ai+A3), A2 = 2(A2 + A3), A3 = 2A3 + A4, 

A4 = -A4 + ^(A5 + A6), A5 = ^(A5-A6). (4) 

In this Letter, we are concerned with the production of a Higgs-boson via 77 scattering. This 
process can proceed through the following two basic and independent mechanisms: 

• Direct scattering of two real photons 77 h, see Fig. [1^; 

• Virtual two-photon fusion in e"^e~ collisions, namely e"'"e~ — > e"'"e~7*7* — > e'^e~ + h (Fig.[T}3). 

Although a tree-level 77/i-coupling is not allowed by the electromagnetic gauge symmetry, this 
interaction is generated at the quantum level through a plethora of radiative corrections, whose 
description in terms of Feynman diagrams is displayed in Figure El The entire set of diagrams 
corresponds to the production process 77 ^ for both CP-even states h = h^,H^, whereas for 
the production of the CP-odd one, 77 AP, only the first line of diagrams is allowed (owing to 
C and CP-invariance) . The relevant cross-section can be written in terms of the partial width of 
the decay /i — > 77 as follows: 

^(77 --h) = —nh^ 77) S{s - M,) (1 + ^1 = 8n ^s-MlY + MlTl ' 

with r/1^2 = il the helicities of the two colliding photons and where we have used a standard 
representation for the (5-distribution. Following the above strategy, and computing the diagrams 
contributing to the amplitude with the help of the computational packages FeynArts, FormCalc 
and LoopTools [26], we may finally arrive at the 77-induced single Higgs boson cross-section. In 
order to obtain more accurate results, a running value for the electromagnetic coupling constant 
aemi^z) = 1/127.9 has been used. 
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Figure 2: Feynman diagrams describing the process 77 — > h°, at the one-loop level, within the 2HDM. In the SM 
case, we have to disregard the diagrams with charged Higgs bosons. 

3 Single Higgs boson production in a 77 collider 

The basic 77 collider operates mainly through the mechanism of Compton (back) scattering of laser 
photons off the linac beams [13]. In a nutshell: a photon from a laser pulse collides with a high 
energy electron (or positron) at a small angle; as a result, the electron recoils and one is left with a 
Compton-scattered photon traveling in the direction of the original incident electron. The efficiency 
of the e"^ ^ J "conversion" will depend on many factors, such as the energy of each of the beams, 
the properties of the laser as well as on several non-linear effects and subleading mechanisms which 
turn out to modulate the overall process. As a matter of fact, the entire procedure comes down 
to furnish a 77 luminosity spectrum, for which several standard parameterizations are available. 
This spectrum is essential to compute the expected number of events (in our case the number of 
produced Higgs bosons of a particular neutral species h = /jO, H^,A^). To this end, the cross-section 
computed in the previous section must be appropriately folded with the normalized (dimensionless) 
photon densities provided by the given parameterizations. The total (unpolarized) 77 cross-section, 
after — ^ 7 "conversion" of the primary linac beam, can finally be engineered from the following 
recipe: 




(6) 
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where the partonic cross-section o"^- r^. is given by 



with s = Ts {s being the CM energy of the primary hnac machine). In the above expression, dC^j/dr 
stands for the (differential) photon luminosity distribution constructed out of the photon densities 

/j/ei,/i/e2 as follows: 



dq^ /■! dx 1 



dr X 1 + 6., 



[/i/ei (^) fj/e2 (t/x) + /j/ei i^) hje^ 'k^' I x)\ ■ (8) 



Functions /j/g^ and /j/e2 (one per beam of given polarization) are taken, in our case, from the 
standard package CompAZ [27]. 

In order to proceed with the numerical analysis, let us first of all examine the structure of the 
partonic cross-section. At fixed s = M|, it can be expressed from Eq. dS]) in a simply manner: 

ct(77 ^ /i) = ^(1 -Fr/i?72)fi(/i ^ 77) . (9) 

If averaged over polarizations, the cross-section is given by the previous result but without the 
factor 1 + ??i because one has to sum over 7/1,772 = ±1 and divide by 4. For polarized photon 

beams of equal polarization (i.e. +-|- or ), the resulting cross-section is a factor of 2 larger. If, 

alternatively, we consider the case of opposite polarizations (H — ), the cross-section vanishes (as 
expected from angular momentum conservation). 

Let us emphasize that the partonic cross-section encodes already the relevant information re- 
garding the dynamical features of the 2IIDM under study. This information is obviously contained 
in the reduced amplitude M of the effective ^177 vertex, and thus also in the branching ratio of 
h — *■ 77- Notice that, in the region where the triple interaction dominates, the amplitude behaves 
roughly as A4(77 "bP) ~ aem Cn+H-ho- Obviously, in this region, the cross-section ([9]) is directly 
sensitive to the H^H^h'^ trilinear Higgs self-interaction (cf. e.g. the first two diagrams of the 
second row in Fig. [2]), given by 



■'H+H-hO 



2Mw sin 



sin(/3 - a) [mIo - 2 M^± 
2M2o -4^A5M^ sin^^H- 



sin 2/3 



(10) 



It follows from this expression that such coupling can be enhanced either at low or high values of 
tan /3, and also through the Higgs boson mass splittings - unlike the MSSM case. In addition, it 
may be heightened through its explicit dependence on the A5 parameter H. Needless to say, in all 
cases these enhancements are strictly harnessed by the various theoretical and phenomenological 
constraints discussed in the previous section. By comparison, the corresponding (approximate) 
behavior of the amplitude for the SM contribution, if we assume that is dominated by the top 
quark Yukawa coupling, reads 7^(77 ^ H) ~ NcQ'laernernf: / {M'jj M\y sin^v^)- This amplitude 
receives in general sizeable contributions from the gauge boson loops that significantly correct 
it. Similarly, in the domains of the 2HDM parameter space where the 3H-coupling (jlOp is not 
overwhelming over the Yukawa and gauge boson couplings (cf. e.g. the first, third and subsequent 



^We note that for A5 = Ae = 2^2 Gf M\o = Ml„/{2 sin^ 6w M^), Eq. (|10|l reduces to the corresponding result 
of Table 1 of Ref. [11], as it should. 
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rows of diagrams in Fig. [5]), our estimate above does not even hold as a crude approximation. 
Therefore, in general it is necessary to come to grips with the full expression for the effective 
coupling gyyh (which is indeed the main object under study, mainly for the cases h = h^,H^). 
Even more useful is to define the exact ratio between the corresponding 2HDM and SM coupling 
strengths at one-loop: 



r = 



\M 



2HDM 



M 



SM 



r{h 77) 
r(/7^77)J 



1/2 



(11) 



Obviously, the 2HDM amplitude will depend on whether the model is of type-I or type-II. In the 
particular region where the trilinear coupling dominates for production, the expectation on 
that ratio is roughly of order r ~ C^+-^-^o Mfj M\y /m^ . Therefore, in the regime of large IA5I, 
for which Cjj+jj-j^o ~ Myy \Xq\, we can foresee a big enhancement with respect to the SM. To be 
sure, in practice we will perform the numerical analysis of the exact expression (jlip and consider 
its behavior in general regions of parameter space. 

For the sake of convenience, let us focus hereafter on the sets of Higgs boson masses displayed 
in Table [T] below. The mass spectrum in Set I of Table [H for instance, allows to enhance the 



2HDM 


Set I 


Set II 


Set III 


Set IV 


Mho 


115 


150 


200 


200 


Mho 


165 


200 


250 


250 


Mao 


100 


110 


290 


340 


Mh± 


105 


105 


300 


350 



Table 1: Higgs mass parameters, in GeV, used throughout the calculation. 

II^II"h° coupling. Its maximum value is roughly attained for sin a = —0.86, tan/? = 1.70 and 
A5 = —25.0. Incidentally, notice that the aforementioned Set I of masses is only suitable for type-I 
2HDM, due to the relatively light value chosen for the charged Higgs boson (Mjji = 105 GeV), 
which is below the (indirect) limit of 295 GeV afflicting the type-II models [20]. Let us recall 
that, experimentally, the current 95% C.L. direct mass limits for general Higgs bosons searches 
are: M^± > 79.3 GeV for the charged Higgs boson, and M^^o > 92.8 GeV, M^o > 93.4 GeV 
(tan/3 > 0.4) for the neutral ones (of course with Mjjo > M^o) [21]. 

In Fig. [3l we perform the numerical analysis of the averaged total cross-section ((t^^_+/i) (s), 
Eq. ([6]), in which the partonic contribution is folded with the effective luminosity function ([8]). 
We display ((T^^^/i)(s) as a function of the center-of-mass (CM) energy y/s of the linac machine 
for the Set I of Higgs masses. In this figure, we explore a region of parameter space where the 
H^H~h'^ coupling dominates for different negative values of A5 (As we will see later on, large 
A5 > values are forbidden by vacuum stability.). Notice that, for sufficiently large [Asl > 10, 
the cross-sections can be considerably high (spanning the range 0.01 — 0.2 pb) at the fiducial 
startup value ^/s = 500 GeV of the ILC, and entailing at this point more than 10^ — 10^ events 
for the given integrated luminosity. The rates, however, decrease fast for smaller values of [Asl, 
the reason being the destructive interference between the charged Higgs, gauge boson and fermion 
loops at low values of [Asl. For example, at the same energy and for A5 = (—2, —5, —8) we obtain 
(a^y^fi) = (2.40,0.16,5.50) fb respectively. Clearly, there is a delicate balance in the low [Asl 
region which amounts to a severe depletion of the overall 2HDM cross-section. By comparison, the 
corresponding result in the SM, namely for the same Higgs mass {Mh = M/jo), reads (c-y-y^//) ~ 11 
fb (hence 5 x 10^ events at that energy and luminosity range), which is quite sizeable. It follows 
that even small departures from this value should be measurable, especially in a high precision 
instrument as a 77 collider. 
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Figure 3: Cross-section {o'-yy^h){s) given by Eq. (|6]) and number of Higgs boson events, as a function of the CM 
energy of the linac, assuming an (integrated) luminosity C — 500 fb~^. We plot the corresponding values for the SM 
and the 2HDM using the Set I of Higgs boson masses, sin a = —0.86, tan/3 = 1.7 and three values for As. 



In Fig.m we test the dependence of the cross-section {a^^^^o) on the Higgs masses. Specifically, 
we plot its evolution in terms of Mj^o (left panel) and Mjj± (right panel). We take the other 
parameters from Set I while keeping sin a = —0.86, tan/? = 1.70, and choose three different values 
for A5: -20, -10 and A5 = Ae = 2V2GfMIo ~ 0.34. Included is also the SM cross-section 
{a-yj^n)- Worth noticing is the fact that, while the 311 coupling is dominant, the cross-section 
does not immediately drop with Mf^o; it actually increases for a while up to a few hundred GeV. 
Moreover, for large values of IA5I, the evolution of {a^^^j^o) presents a notorious "spike-shaped" 
enhancement near M^o ~ 2 • Mjji ~ 210 GeV, which is brought about by the threshold effect of 
two real charged Higgs bosons in the loop. The corresponding effect for a pair of real vector bosons 
at MijO ~ 2 • My/± ~ 160 GeV is also barely visible therein. 

The sharp suppression dip standing out on the right panel of Fig. d] deserves also a few words. 
Recall that, in the regime in which the influence of the trilinear contribution (H^H^h*^ ) greatly 
"waxes", the charged-Higgs mediated correction holds absolute sway over the loop-induced cou- 
pling g^yho- However, as soon as we raise the charged Higgs mass, the positive influence of the 
trilinear Higgs boson interaction rapidly wanes and it cancels more and more against the (negative) 
effects from the gauge boson and fermion loops. This gives rise to the aforementioned destruc- 
tive interference. Eventually, a particular value of M-^± is reached where the two sorts of effects 
virtually annihilate each other (right at the vertex of the dip in the figure). Beyond this point, 
one rapidly reaches a regime where only the gauge boson and Yukawa coupling (negative) effects 
remain. Most of this region is actually excluded by the constraints. The destructive interference 
pattern described here is only possible when the set of Higgs boson masses is relatively light, as in 
the case of Set I under consideration, otherwise the trilinear effects could not be competitive. 

The behavior of the effective coupling g^^ho in the 2HDM can be better assessed in terms of 
the ratio r defined in (jlip . Its dependence on the Higgs mass spectrum is sampled in Table [2] using 
the parameter setups indicated in Table [TJ Set II, for instance, contains a heavier neutral CP-even 
Higgs sector which results in a sizeable value of r of 3.75. It means that, in this case, the effective 
strength of the 77/1*^ vertex almost quadruples that of the SM ('yyH). This is quite remarkable. 
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Figure 4: Cross-section {cr^^_;jo)(s) as a function of the light CP-even Higgs mass (left panel) and the charged 
Higgs mass (right panel) at fixed ^/s = 500 GeV. Remaining mass parameters from Set I, sin a = —0.86, tan/3 = 1.70 
and three values for A5: —20, —10 and A5 = Ae- The SM cross-section (cr^^^^o) is also included (in the right panel, 
it almost coincides with the A5 = Ag case) . The shaded area is excluded by the constraints. 



On the other hand, Sets III and IV are characterized by heavier charged Higgs bosons and at the 
same time by a heavier CP-odd Higgs bosons (so as to elude the 5p bounds). Unsurprisingly, the 
ratio r falls in this case to within values below 1, i.e. close to the SM from below. As it should be 
expected, the larger the charged Higgs mass is, the less efficient is the enhancement capabilities 
associated to the 3H self-interactions. Incidentally, let us notice that Sets III and IV of Higgs- 
boson masses are intended to describe type-II 2HDM. Does this mean that for type-II models 
(those closer to the MSSM Higgs sector) there is no hope to hint at non-SM Higgs boson physics 
with 77 collisions? Not necessarily so, as there is a tail of subleading one-loop effects triggered by 
the non-SM Yukawa couplings of the 2HDM in combination with the gauge bosons; in particular, 
we have already detected it in Fig. [Ha for Set I (although within the excluded region). But, in 
general, this tail is available and lies well within the allowed region for heavier sets of Higss boson 
masses, such as Sets III and IV. More on it below. 

A brief comparison with the existing calculations of single Higgs boson production within the 
MSSM is in order [17, 18]. For example, in Ref. [18] the ratio between the decay widths of the 
MSSM Higgs boson and the SM Higgs boson H into 77 is computed at one-loop. This ratio 
corresponds to the square of r defined in (jlip. They find that, in the most favorable regimes 
(viz. large mass-splitting and large mixing angle in the top squark sector), it can lead to values 
of r up to \/2 ~ 1.4. It follows that the most optimistic MSSM expectations on r are markedly 
below the maximum enhancement capabilities of the 2HDM (r ~ 4). The reason for this is 
clear and it was already advanced in the introduction - namely, the 3H self-couplings within the 
MSSM are restricted to be gauge-like and cannot source the potentially large effects that we have 
identified within the general 2HDM. Therefore, what we have called the "tail of subleading effects" 
in the 2HDM case is actually one of the main sources of the MSSM effects, the other being the 
supersymmetric Yukawa couplings of the Higgs bosons with the squarks. 

Figure [5] presents the corresponding contour lines for the ratio (jlip in the (AsjSina) plane. 
The cross symbol indicated on the left at the bottom denotes the point in this plane where such 
ratio attains the maximum value permitted by all the constraints. Moreover, we show the regions 



9 





Set I 


Set II 


Set III 


Set IV 


r 


3.98 


3.75 


0.98 


0.98 


tan/3 


1.7 


1.7 


1.0 


1.0 


sin a 


-0.86 


-0.86 


-0.82 


-0.82 


As 


-25 


-25 









Table 2: Maximum value of the ratio r, Eq. (Illf) . in the case of /i° and for the different mass sets quoted in Table[l] 
together with the configuration of tan/3, sin a, A5 for which these optimal values are attained. 




-1 -0.5 0.5 1 

sin a 



Figure 5: Contour hues for the ratio r — g^-yh" /djjH (HJ) in the (AsjSina) plane, for the Set I and tan/3 = 1.70. 
The regions excluded by vacuum stability (upper half plane A5 > 0) and the unitarity of the two kinds of Higgs boson 
self-couplings, namely the quartic (large circular domains) and the trilinear (small circular domains), are separately 
shown. The cross symbol on the left (at the bottom) denotes the point with maximum allowed r. 



excluded by vacuum stability and by unitarity of the trilinear and quartic Higgs boson couplings. 
Notice that the sign A5 > is mostly forbidden by vacuum stability, which explains why we have 
presented the numerical analysis of the previous figures only for A5 < 0. In any case, we see that 
there is a sizable region left where the effective 2HDM coupling g^^ho is significantly larger (in 
fact, a few times larger) than the SM coupling g^yyH- 

Finally, the available domains across the (A5,sina) plane wherein one can obtain enhanced 
values of the ratio r with respect to the SM case is explored systematically in Fig. [6] for different 
values of tan/3 and at the fiducial startup energy y/s = 500 GeV of the ILC. Specifically, in this 
figure we compute those regions in which the predicted value of the ratio r exceeds the correspon- 
ding SM value (r = 1) by 10% at least while still being compliant with the full set of constraints; 
equivalently, regions where the cross-section is augmented by 20% or more, thus inducing an excess 
of about one thousand events above the SM prediction (within the given luminosity segment). As 
we can see, the projected domains are sizeable. Interestingly enough, even tiny ±1% departures 
of r from 1 would already be significant, as they would amount to 100 events up or down the 
SM expectation. If we would adopt this softer criterion, the allowed domains in Fig. [6] would be 
substantially augmented. 

Particularly stringent is the impact of the unitarity restrictions, which translates into a pro- 
gressive shrinking of the relevant region as we raise the value of tan /3. For completeness, we have 
also addressed the single production of a heavy CP-even and a CP-odd Higgs boson. Their re- 
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Figure 6: Regions in the (As, sin a) plane (allowed by all the constraints) where the ratio is r > 1.1, i.e. when 
Qj^h is at least 10% bigger than g-y^H for tan/3 = 1.3, 1.7 and 2.1 



spective maximum cross-sections (at ^/s = 500 GeV) for the specific case of Set I read as follows: 
{cr^^^H')) = 0.31 pb and ((T^^__^^o) = 1.9 fb. The former is of the order of {a^^^f^a) itself, whereas 
the latter is significantly smaller, the depletion being caused by the absence of trilinear couplings 
of with the charged Higgs bosons. Furthermore, we wish to emphasize the existence of a re- 
gion of parameter space where {a^^^uo) and {a^^^f^o) can both be simultaneously sizeable. That 
region (e.g. for Set I) is just the allowed domain in the down-right corner in Fig.[5l There we find 
— 0-2 pb and (c7^^_>/jo) ~ 0.05 pb, entailing some 10^ — 10^ events within the standard 
luminosity range. This possibility can be very relevant, as it could be responsible for a potentially 
distinctive 2HDM signature which is unmatched in the MSSM. 

We have also tested the enhancement potential of the 2HDM in the domains of parameter space 
where the 3H-coupling (jlOp is not relevant and where the bulk of the contribution concentrates on 
the Higgs-fermion Yukawa couplings in combination with the gauge bosons. In the case of generic 
type-II models, there is no enhancement at large tan (3 (unlike the supersymmetric case) inasmuch 
as tan (3 is severely constrained by unitarity. In general, due to the destructive interference between 
the diagrams dominated by the SH-coupling and the rest (fermion and gauge boson loops), the 
enhancement capabilities of the Yukawa sector become overshadowed. As a consequence, in such 
region we meet the following situation: 1) There are still non-negligible domains in the 2IIDM 
parameter space where the cross-section departs remarkably from its SM counterpart. However, in 
most cases the departure entails a significant (e.g. 10%) reduction of the cross-section with respect 
to the SM; 2) type-I and type-II models become essentially indistinguishable in that domain. This 
is a reflect of the fact that the Higgs-top quark coupling (which has the same form in either type-I 
and type-II models) drives the leading contribution in the Yukawa sector, whereas the gauge boson 
contribution is common in both types of models. Therefore, spotting a tail of non-SM effects in 
this region could not distinguish the type of 2IIDM. Still, the missing number of events could 
certainly hint at the existence of a smoking gun triggered by physics beyond the SM. 

Let us close this section by briefly mentioning that we have also surveyed the impact of another, 
more restrictive, set of unitarity constraints [22]. We have found that, in the most optimistic 
scenario for h*^ production in 77 collisions, the II'''H~h^ coupling lies roughly a factor 2 — 3 below 
the largest value it can take under the current set of constraints ([1]) and ([2]). Correspondingly, 
the enhancement with respect to the SM would be a factor 5 — 10 times milder, thus dwarfing the 
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relevance of the main effects in some regions of the parameter space. We point out, though, that 
the unitarity restrictions proposed in [22] are not fully coincident with those considered in [23], 
and in this sense there is still some controversy in the literature on this issue. 



4 Single Higgs boson production through 77 fusion 

The interest on virtual photon- fusion ( "two-photon processes" ) certainly has a long and widespread 
history in Particle Physics. For instance, a prospect for the measurement of the pion lifetime from 

e+e" ^ 7*7* ^ e+e" + 7r° (12) 

was first discussed by F. Low almost half a century ago [28]. Studies of these processes, as well as 
detailed considerations on two-photon production of muon pairs and multi-pion final states, were 
carried out subsequently in the seventies [29]. Remarkably enough, long after the first pioneering 
studies appeared, two-photon processes are still an active and fruitful field of investigation, in 
particular for Higgs boson production. In actual fact, single Higgs boson production is, in a sense, 
the modern counterpart of Low's "single-meson" production from quasi-real two-photon collisions. 
Already in the early eighties, this Higgs production channel was first studied in the literature 
within the context of the old e^e~ (pre-LEP) colliders by J. A. Grifols and R. Pascual [15]. 

To be sure, the traditional two-photon processes are the forerunner of the future 77 colliders 
considered in the previous section. These colliders will probably concentrate most of the interest 
around linac physics in the future and may greatly supersede the former in all practical searches 
for new physics. It is instructive to see once more why, specially in regard to the "Higgs issue", a 
most sensitive matter these days. To this end, we compute here the cross-section for the processes 

e+e" ^ 7*7* ^ e+e" + h {h = h°, H°, A^) , (13) 

and compare with the results for real 77 collisions studied in the previous section. In practice, 
we concentrate on the lightest CP-even state. Noteworthy is the fact that the cross-section for 
the virtual 7*7*-fusion processes, in contradistinction to the real 77 collisions, grows with the CM 
energy up to very high values of ^/s. Generically, the behavior of e'^e" ^ 7* 7* Y + e~^e~ in 
the asymptotic energy regime goes as ~ (a^/A/P) ln^{s/ml) ln"(s/M2), where M is the threshold 
mass of the produced final state Y, and the number n > 1 depends on the high energy behavior 
of 17(77 The logarithmic growth simply tracks the dynamical feature by which, for these 

processes, the virtual photons 7* can be quasi-real and hence have their momenta well-below the 
CM energy of the process, which may satisfy s ^ My - the rest of the energy being carried away 
by the concomitant lepton final states. Recently, the pairwise production of Higgs-bosons via weak 
gauge-boson fusion mechanism, e+e~ V*V* hh + X {V = W"^, Z\ h = h°, H°, A°, H='=) was 
analyzed in Ref. [12] at the leading order 0{a'^^). The process turned out to be instrumental 
for probing the 3H self-interactions. A complementary strategy along the same lines is offered 
by the 77-fusion process (fT3|) . in which the loop- induced 77/1 vertex at order 0{a^^aew) can be 
dominated by the 3H coupling. 

The technical complications associated to two-photon processes of the kind ()12p -()13 p were 
tackled long ago in the literature (cf. [29] for a classical review). The physics is also well understood, 
it boils down to the well-known equivalent-photon or Weizsdcker- Williams approximation, by which 
the virtual photon emitted from the scattered electron appears near the mass shell. This feature 
allows to essentially trade the process of electro-production for a photo-production one with an 
appropriate photon spectrum. The so-called photon content of a given electron can be explicitly 
factorized and the production cross-section can be approximated in the following way: 

t V / .S I 

cT(e+e- ^ e+e-X) 



) J f{r)(J^^/^x{Ts)dT, (14) 
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Figure 7: Evolution of the two-photon cross-section cr(e"'"e~ e^e~h'') as a function of the CM energy and the 
corresponding number of events for a total luminosity C = 500 fb^^. We include the 2HDM and SM curves using 
the Set I mass parameters, sin a — —0.86, tan/3 — 1.7 and four values of A5. 



where tq = M'^/s and /(r) = (1/r) [{2 + log(l/r) - 2 (1 - r)(3 + r)] . 

Figure [7] presents the logarithmic evolution of the production cross-section cr(e"'"e~ — > e"'"e~h'') 
as a function of the CM energy. We have used the mass Set I and the optimal parameters quoted 
in Table [2j Once the production threshold has been surpassed, the cross-section suddenly increases 
up to a value of about 10^^ pb and becomes persistently sustained, with only a mild (logarithmic) 
evolution. Over this approximate plateau, it amounts to more than 5000 events for an integrated 
luminosity C = 500 fb~^, therefore resulting in an enhancement of one order of magnitude with 
respect to the SM. Even more remarkable is the fact that, despite the fairly large value of the 
two-photon cross-section, it is still a factor of 10 (at least) smaller than that of the corresponding 
77 real scattering processes, for the same values of the 2HDM parameters (cf. Fig [3]). This should 
explain vividly and convincingly the outstanding superiority of the future 77 colliders versus the 
ancient two-photon processes. 

5 Conclusions 

We have devoted this work to analyze the production of a single neutral Higgs boson, h = 
h^, H^, A^, within the two-Higgs-doublet model (2HDM) through the following complementary 
mechanisms: i) direct scattering of real photons in a 77 collider; and ii) fusion of a virtual photon 
pair in a conventional two-photon process, e"'"e~ — > e^e~7*7* — > e"''e~ -|- h. Both mechanisms 
are direct handles on the effective 77/1 interaction, g-y-yh- This coupling is a pure quantum effect 
generated by a plethora of radiative corrections involving charged Higgs bosons, quarks and gauge 
bosons. Among the interactions in the loops, we have the 3H self-interactions, most remarkably 
H'''H~h'' in the case of production - which proves to be utterly dominant in certain regions of 
the 2HDM parameter space. We have systematically swept this space and identified those configu- 
rations for which the departure from the SM prediction is most remarkable, and we have done this 
in full compliance with the rigorous constraints dictated by perturbativity, unitarity and vacuum 
stability bounds, as well as by the EW precision data. The result is that, in the most favorable 
scenarios, the physical cross-section (i.e. the one convoluted with the backscattered luminosity 
function) can typically reach the level of (cr-yy^h) ~ 0.1 — Ipb. In other words, (a-y-y-^h) may 
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rocket to values 10 — 100 bigger than the expected SM yield (a^^^n) (for similar values of the 
Higgs mass), which by itself should already be perfectly measurable: (a^^^n) ~ 0.01 pb = 10 
fb. Such notorious enhancement can be traced back to the behavior of the H"'"H~h'' coupling, 
despite it is highly restrained by the overall constraints. By sticking to moderate tan/3 > 1, it is 
possible to licitly increase the rates by choosing relatively large (negative) values of the parameter 
A5. Moreover, in order to optimize this mechanism, the charged Higgs boson should be relatively 
light (say, below 300 GeV) so that the associated quantum corrections are not severely hampered 
by the decoupling effects. It means that the ^^^/j-enhancements that we have encountered apply 
only for type-I 2HDM, because for this kind of models the mass of the charged Higgs boson is not 
constrained by B(b sj). A similar conclusion ensues for the case of double Higgs production in 
77 collisions, but with lower cross-sections [14]. 

In the above conditions, the expected number of single Higgs boson events emerging from direct 
77 /i scattering, within the typical energy range of the ILC (500 — 1000 GeV), is of the order of 
10^ per 500 fb"^ of integrated luminosity. Compared to the production rate of single Higgs bosons 
through the traditional virtual photon-pair fusion, e+e~ e"'"e~7*7* — > e+e" -|- h (which performs 
at the level of 0.01 pb, at most, for ^/s > 500 GeV), the real 77-collision mechanism is at least one 
order of magnitude more efficient. 

On the experimental side, the prospects for Higgs boson detection in a 77-collider are deemed 
to be excellent. To start with, let us stress that the single Higgs-boson final state is to be produced 
essentially at rest. Therefore, for < 2My < 180 GeF, the corresponding signatures should 
mostly be in the form of back-to-back, highly energetic, quark jets (66, cc). For > 2Mv, instead, 
signatures with two or four charged leptons in the final state (from ^ £^ + missing energy 
and, specially, from Z i~^£~) should be really pristine. Furthermore, we have seen that, in some 
cases, the two channels 77 — *• hP and 77 — are simultaneously accessible and with similar 
rates. Needless to say, this could result in a double distinctive signature of new physics. 

With enough statistics on these events, and upon analyzing the invariant mass distribution of 
the resulting jets of quarks and leptons, the measurement of some 2HDM Higgs boson(s) mass(es) 
should be attainable with fairly good accuracy, together with a precise determination of the effective 
Qjjh couplings (typically for hP and/or H^). If their strengths would happen to be vastly dominated 
by the triple Higgs boson self-interactions, the signature of non-standard Higgs boson physics would 
be crystal-clear, leading us to suspect it to be rooted in some generic type-I 2HDM. However, should 
we meet the juncture g^^h ^ d'yyH or g^^^ > g^^yH, the underlying quantum effects would be largely 
insensitive to the type of model and a more detailed comparative study with the MSSM would 
be mandatory [30]. Even then, tiny deviations could hint at new physics. In this regard, it is 
important to emphasize that, given the high precision nature of a 77 collider, gathering a small 
5 — 10% effect (positive or negative) should be sufficient to point at a smoking gun. 
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